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SU^WARY 


The linearized equations describing the propagation of sound 
in variable area ducts containing flow are shown to be singular 
when the duct mean flow is sonic. The singularity is removed when 
previously ignored nonlinear terms are retained. 

The results of a numerical study, for the case of plane waves 
propagating in a one-dimensional converging-diverging duct, show 
that the sound field is adequately described by the linearized 
equationsonly when the axial mean flow Mach number at the duct 
throat M^j^<0.6. For numerical results showed that 

acoustic energy flux was not conserved. An attempt was made to 
extend t.'.e study to Include the nonlinear behavior of the sound 
field. Meaningful results were not obtained due, primarily, to 
numerical difficulties. 



DEFINITION OF SYMBOLS 


Symbol 

A(x) 

c(x) 


G 

i 

2L 

M 

X 

m' 

P.P' 

S 

R 


r 

T 

t 

u' 

V 

X 

a 

B 

Y.fi 

e 

c 

X 


Definition 

duct cross-sectional area at location x 
speed of sound at location x 

non-dimensionalized acoustic density and velocity 
functions respectively defined by Eqn. (40) 

imaginary number fT 

sound source wave number (^w^/c^) 

duct converging-diverging length 

duct mean Mach number (=V^/c) 

acoustic Mach number (=u^'/c) 

mean and acoustic pressure respectively 

duct vector velocity (=V+u' ) 

sound energy reflection coefficient defined by 
Eq. (69) 

cylindrical radial coordinate 

sound energy transmission coefficient defined 
by Eq. (69) 

time 

acoustic velocity 
duct mean velocity 
duct axial coordinate 
see Eq. (24); also Eq. (62) 
see Eq. (25), also (62) 
defined by Eq. (62) 

a cylindrical azimuthal coordinate 

I d(^ 

G(k) 

G(x=o)x^G(x) 

ii 



Symbol 


Definition 


\> 


\h/% 


p,p' 


mean and acoustic density 


£ 


sound radian frequency 
defined by Eq. (62) 


V 

Subscripts 


phase angle 

column matrix defined by Eqn. (21) 


r,9,x 

denotes radial, azimuthal, or axial cylindrical 
coordinates 

o,in 

denotes 

duct section at x = 0 

2 

denotes 

nondiraensionalized duct section at x = 2 

r,t 

denotes sound reflection and transmission 
respectively 

p,u 

denotes 

acoustic pressure and velocity respectively 

s 

denotes 

sound source 

th 

denotes 

duct throat 

T 

denotes 

total or stagnation values 

(J 

denotes 

vector quantity 

R.I 

Superscripts 

denotes 

real and imaginary parts respectively 

(“) 

denotes 

matrix quantity 

n 

denotes 

mean value 


111 



1.0 INTRODUCTION 


Jet engine inlet whine is a particularly troublesome and 
persistent source of community noise pollution. A promising new 
technique to reduce or eliminate this noise is the use of variable 
area inlets to choke the flow. The idea here is that the choked 
flow acts as an effective barrier to sound transmission through the 
inlet. Recent experimental sonic inlet noise suppression studies^*^ 
demonstrate the feasibility of this concept. The success of these 
studies have prompted a number of theoretical studies dealing with 
sound propagation through variable area ducts containing flow**®. 

All of these studies, however, share a common deficiency - they 
assume that the usual linearized convected wave equations adequately 
describe the bahavior of sound waves in an accelerating flow. 

The principal objective of this study is to show that the 
linearized sound propagation equations are not valid for sonic or 
near sonic mean flows. The linearized sound propagation equations 
are singular when the mean flow is sonic. The singularity 
is removed only when nonlinear terms previously ignored are 
retained. A second objective is to understand the details of how 
the sonic inlet reduces the transmission of sound through the inlet. 
As we understand it, there are two effects that act to reduce the 
sound transmission. One is related to the fluid mean velocity at 
the throat (called convection) and the other to flow inhomogeneities 
(i.e., mean velocity and density gradients). 


The reduction of inlet sound transmission by convection is 
explained as follows. The rate at which sound energy propagates 
through a duct containing inhomogeneous media is called the group 
velocity V . For simplicity consider a highly idealized sonic inlet 
consisting^of a constant area duct containing a uniform flow. For 
this case, the group velocity simplifies to 


V = c-V 
g 

where c is the local speed of sound and V is the duct mean flow. 

It is clear that the rate at which acoustic energy is transmitted 
out of the inlet decreases with increasing speed. When V=c, 
there will be no transmission of sound out the inlet; all the inter- 
nally generated sound will be reflected back into the sound generat- 
ing interior region of the duct. A physical explanation of the 
breakdown of the linearized sound propagation equation when the duct 
mean flow is sonic can be made in terms of group velocity. 

The reduction of the rate at which the sound energy is transmitted 
out the duct inlet can be thought of as a piling-up of the sound 
at the throat. When the duct flow is choked, all of the sound piles- 
up; thus the sound pressure and density increase at the throat 
violating the original linearized assumptions of small disturbances. 



It follows that under these circumstances nonlinear terms, previously 
ignored, must be retained. 

When the flow contains inhoraogeneities in the form of mean 
velocity and density gradients, the situation becomes very complicated. 
It is no longer easy to describe analytically the group velocity. 

The distinction is made here between gradients that cause refraction 
and gradients that cause reflection (or scattering). Gradients 
perpendicular to the wave vector refract or bend the wave and grad- 
ients that are in the same direction as the wave vector reflect or 
scatter the wave. In this study, only scattering by mean axial grad- 
ients will be considered. 

The report is organized as follows. The behavior of sound 
waves propagating in a converging-diverging choked duct are described 
in Section 2, The Mach number range within which the linearized sound 
propagation equations are valid is explored numerically in Section 3. 
The report closes with a summary of the results of the study and 
describes aspects of the behavior of sound fields in accelerating 
duct flows requiring further study. 

2.0 GOVERNING SOUND PROPAGATION EQUATIONS 

The equations describing the propagation of sound in a variable 
area duct containing a mean flow will be derived in its most general 
form. The objective here is to show that the linearized sound propa- 
gation equations are singular when the duct mean flow is sonic. The 
(variable area) duct mean flow V is assumed to contain cylindrical 
components (V ,Vg,V ) where the subscripts (r,6,x) represent the 
radial, azimufhal, ind axial directions respectively as defined in 
Figure 1. In general, V is a function of all three components. The 
flow is assumed to be inviscid and nonheat- conducting so that the fluid 
pressure and density are adiabatically related. Further, the mean 
flow is assumed to be independent of the sound field. In the follow- 
ing derivation the quantities (p,p,c) represent mean fluid pressure, 
density and sound speed and (p',p’,u') represent acoustic pressure, 
density, and velocity respectively. Vector quantities will be denoted 
by a single horizontal bar located below the symbol ( ) and matrix 
quantities by a double horizontal bar located above tKe symbol (“) . 

The equations describing the conservation of mass and momentum 
are respectively 


iL 


+ V-(p9)*0 


and, in component form, 


^ bOr . 5®. IPr 
^ de 



( 1 ) 


( 2 ) 


2 



+ Qr + -Qo ^^9 4 . QrQe ^ Qj( ^P = q 

p(j^^ + Qr AQx + ^ ^Qx -f- Qxj^Qx^ + ^P r=n 

^V dt ar r 5^ b^JhK 


where 


Q =^Qr, Qe, Qx) ~ Y 


is the total flow velocity. In Eqns. (3) -(5), the variations of 
pressure was replaced by the adiabatic condition that 


%= 


where c is the fluid local sound speed. Now assume that the flow 
field can be written as the sum of the mean and fluctuating quantities 
so that 


p(r,®,y,tj =p(r,e^x)+ p' (r,6,x,-t) 
p(r,©,x,t) = p(r,6,x)+ p'(r,a,x,-t) 




It is further assumed that the time-averages 


/ / 


P =P 


= u' = 0 


3 



1 f 

where ( ) ~ f Jo ( )dt and T is the sound wave period. Substituting 
Eqns. (7-9) *and retaining all terms yields 

> 

(v-y)p'+(fy)(v-w')+(H-v)(f+p)+(v-v p' = 0 


dt 


(11) 


and 


iHC +(u'v)(Vr+Mr)+(y.VjM'r--L(2Ve+ Ue)ue 



( 12 ) 


+(w'- v)(Vet tte) + (y. V)t< 0 t(Vr + Wr)l 

r r 





= 0 


+(«'. v)(Vx+Mx)+(y. v)m,' + c’(i- 



(13) 


(14) 


where the operator V equals 


V=/— J- 

\(3r ^ r ^ d% / 


(15) 


The notation used to express Eqns. (11-14) is considerably simplified 
using matrix notation. Letting 


4 



1 

hr 

-■& 

<)Vr 

<)x 

v^; ar 

X^CrVe) 


^Ve 

hK 

r\^/ 


1 


-(LtiL 


^ he 

hx 

'^e' a* 

i£ 

dr 

_L ^ 

iL 

hx 

V- (y+u) 

(Vp+i^r) 

0 

0 


0 

(vv + Mp) 

0 

0 

0 

0 

(Vp+Wp) 

0 

P 

0 

0 

(Vp+M*^ 


0 

0 

0 

0 

('^e+“$ 

0 

f('i) 

0 

0 

(X9 + ^©) 

0 

0 

P 

0 

(Vq+Wq^ 


( 16 ) 


( 17 ) 


( 18 ) 


5 



and 



0 0 



0 


0 

0 


0 

0 



e 



09 ) 


Substituting the matrices defined by Eqns. (16-19) into Eqns . (11-14), 
the governing equations become 


= A 

hi 


+ M, + 


.=* A 

hr 




( 20 ) 


where I is the identity matrix and $ is the column matrix representing 
the acoust i c field de fined below as 




( 21 ) 


Equation (20) may be rewritten as 



= v-l 


-(M4) 


hi hi r ^ 


( 22 ) 


where the inverse matrix (NU)"’ 


can be shown to equal 


6 




Here the quantities a and B are defined as 


«l=Vx+W,' 


C24) 


P = '-^ (25) 

Ke are now in a position to demonstrate a principal conclusion 
of this report - namely, that the well known linearized wave propa^’a- 
tion equations that are often used to model the behavior of sound 
waves in choked duct^ become singular when the duct mean flow is 
sonic. The matrix (Mn)*^ in Eqn. (24) contains the terms [o^ (a*-c^B) ]"* , 
Substituting from Eqns. (24) ^ (2S) for the quantities a and 3 yields 



where M = V /c and m' = u ’/c. Now recall that in the derivation 
of Eqp.^(2i), all terms were retained including the nonlinear terms. 
If the nonlinear terms were ignored, there the above expression 
simplifies to 

which is clearly singular at = 1 . It is also clear that the 
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nonlinear terms must be retained at = 1 to remove the singularity. 
This conclusion is quite general because the flow fielC within the 
duct is fully three-dimensional and unsteady. 


3.0 NUMERICAL EXAMPLE 

A numerical study has been undertaken to map out, for a particu- 
lar duct geometry, the Mach number range for which the linearize:^ 
sound propagation equations are valid. To simplify the analysis, 
only plane waves propagating in a one-dimensional duct are considered. 
The duct geometry, shown in Figure 2, consists of two infinitely long 
ducts of equal area connected by a contraction-expansion (constriction) 
region. The use of infinitely long ducts avoids the complexities 
associated with end reflections. Mean flow is introduced into the 
duct as shown and is accelerated to sonic or near sonic speeds in the 
throat. Sound propagating in the direction of the mean flow ' 

Fig. 2a) is called the downstream case and sound propagating against 
the mean flow (see Fig. 2b) is the upstream case. 

The duct mean flow properties are assumed to be steady, one- 
dimensional and isentropic. They are described by known functions 
of the duct axial coordinate x, the origin of which is chosen as 
shown in Fig. 2. The constriction region is specified to be symmetri- 
cal about the throat, the area change occurring over the distance 2L. 
The duct cross-sectional area is described by the expression 


a(x ) 

A(oj 







(26) 


where v = A(L)/A(0) is the ratio of the du' t constriction throat area 
(at X = L) to inlet area (at x * 0) . 

3. 1 Derivation of Governing Equations 

The governing linearized sound propagation equations describing 
the behavior of the sound field in the duct are derived below starting 
from the basic mass and momentum conservation equations. Although 
they could be derived starting from the more general form given by 
Eqn. (22), the highly simplifying assumptions of a one-dimensional 
duct and plane sound waves motivates the rederivation of the governing 
equations. 
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The derivation starts from the equations describing the con- 
servation of mass and momentum in a one-dimensional variable area 
duct 


i£. 




(27) 


and 






(28) 


where p, p and Q ere the one-dimensional counterparts of the three- 
dimensional flow field described in Section 2. Assuming that the 
flow field can be described by a steady state and a fluctuating part, 
the equations describing the conservation of mass and momentum simplify 
to 






and 


V^p'+ sV^' + IE-' =0 

r ^ r dx dx r r dx 


(30) 


Implicit in this derivation is that the steady state terms satisfy 
independently of the acoustic terms the time-average conservation 
of mass and momentum. Equations (29) and (30) are further simplified 
by incorporating the following 




(31) 


( 32 ) 


9 



and 





1 = _ 0-M^) JM 

A dx GM dx 


( 33 ) 


(34) 


where the subscript T denotes fluid stagnation or total values. Equa- 
tion (31) represents the adiabatic relationship between the sound 
pressure and density while Equations (32), (33) and (34) follow from 
one-dimensional steady state, isentropic flow. Substituting Eqns. 
(31-34) into Eqns. (29) and (30) yields 


it' 

a 


-K 


Pt 


Pt dM 1 

U 4. 

1 

O 

1 



eC4-i) 


M dx 


1 

<3 

1 

dx 

G% dx 


(35) 


and 




o 

1 

dvL 
- ^ 

Ct 

dM 

1 

VL 4 . 


[G'/2 

dx 

6^- 

dx 


.fT 


(—) 

dx 




Equations (35) and (36) are non-dimensionalized by scaling 
p’, u', X and t as follows. Let 


e=e.f 




X = K,"' 1C* i t = 


-I .It 

CA>* t 


(37) 


where Pg and Ug represent the magnitude of the known acoustic source 
density and velocity introduced into the duct far upstream of the 
constriction. Kg is the acoustic wave number (Kg = Wg/c’o) where Ug is 
the sound radian frequency and Cq is the local speed of sound at x = 0. 
Substituting Eqn. (37) into Eqns. (35) and (36) yields 
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( 38 ) 



and 



+ 


\ G / G d 



= 0 


(39) 


To solve Eqns. (38) and (39), the following form of the solutions are 
assumed, 



(40) 


Substituting Eqn. (40) into Eqns. (38) and (39) yields 


V / 'Q/ Q dx* ^ ( 3 M dx'* 


and 


6« \ r-i d-f * 

dx* 




= 0 


(41) 



-f— ) 

\ W-i', 


3 


(42) 
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The non 'dimensional functions f and g (the ( )* notation is deleted 
for convenien-e) are complex and are written 


^(x) = -?,(*)+•■+,(*) 


( 45 ) 


and 


g(x) = 3*(x)+i 3 j(x) 


(44) 


Substituting Eqns. (43) and (44) into Eqns . (41) and (42) yields 
for the derivatives 


^ - 
da* 

dx 


-ocM^(l-hA) 

-MV** 

-2oi X 


3« 

d9i 

dx 

1 



(2zX.'\ 

-2othA5i 



dx 

|-M‘ 

0 


_l/ 

oc(l + M^) “MA 



dx 



0 

MA^ 0t(l-4-M*) 




where X ~ Q ** I (x) , I M ^(o) , and <jc= a x ■ 

Here again, Eqn. (45) is singular at M = 1. 


3.2 Boundary Conditions 

The boundary conditions are particularly simple because the duct 
cross-section and mean flow are uniform for |x| > L. The boundary 

conditions are different for the cases of upstream and downstream sound 
propagation. 

3.2.1 Downstream Propagation 

At X = 0 (see Fig. 2a), both incident and reflected waves exist. 
The incident wave is assumed known and generated far upstream of the 
constriction. The reflected wave, generated by the mean flow inhomo- 
geneities in the duct constriction is unknown. At x = 2 , only 
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transmitted waves propagating to the right (i.e., right running waves) 
are assumed. 


For values of x < 0, both the duct cross-sectional area and the 
mean flow velocity are constant. Thus the sound propagation equations 
reduce to the well-known (dimensional) solutions 







(46) 


and 




< e 





VLr e 




(47) 


Here pg and us are the assumed kno%m source strengths of the incident 
acoustic density p' and velocity u' respectively, p'j. and Ur are the 
unknown strengths of the reflected sound wave, and is the unknown 
phase of the reflected wave relative to the incident wave. 


Equations (46) and (47) may be suitably non-dimensionalized 
using Eqn. (37) yielding for the (non-dimensional) functions f and g. 


and 


where 


= 


IX 

l-h«V4 


r. 


iX 

l-M 





(48) 


(49) 


(50) 


Equation (50) follows by substituting the solution given by Eqns. (48) 
and (49) into Eqn. (29) (the linearized continuity Eqn.). Thus the 
boundary conditions at x = 0 are 


-P(o) - l+FoC " = [^l + r« cosH',] +»[Fo sinV^ 


(51) 
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and 


3(°) = • ~ r# e -To CosH*,j - i j^r* s«n<P,j 


(52) 


For values of x 2L, (see Fig. 2b), only transmitted (right 
running waves) are permitted. They are described by the solution 


(>*0'.+) = ft « 




( 53 ) 


and 


e 




(54) 


where p'^, u'^ are constants (it can be shown that u'^ by 
substituting Eqn. (53) and (54) into Eqn. (29)) and is a unknown 
phase shift due to nonhomogeneities generated by the duct constriction. 
Non-dimensionalizing Eqns. (53) and (54) using Eqn. (37) yields the 
boundary condition at x = 2L (or x* = 2 K 5 L) 

-P(2KsL) = g(zKjL) = cosf^ + I r^sinf^ (55) 


where 



2KtL 
I + M(o) 


(56) 


3.2.2 Upstream Propagation 


The boundary conditions for upstream sound propagation are very 
similar to the downstream case. Here, the sound approaches the con- 
striction as shown in Fig. 2b. At x = 2L, it is straight-forward to 
show that the incident and reflected wave may be written 




2K,L 

l-M. 


+ 




iSh.) 


( 57 ) 


V. 


Incident Wave 


Rcplccted Wave 
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and 


j(2K,L)=-e I -Mo + Fa e ^ l-M* '' 


^ v_ 


Incident Wave 


Reflected Wave 


at X = 0, the transmitted wave is 

-f (o) = - . = (F„ COS to) + I (Fo sin to) 


(58) 


(59) 


Thus there are four equations gn, gj, fj and four unknowns 

( a » • : >';'o , -r ) , the solution or which completely describes the 
problem. 


5.5 Nume r ical Integration Scheme 


The solution to the four simultaneous differential equations 
described by Eqn. (45) requires that four (4) constants (of integra- 
tion) be specified. The four unknowns (To,^ , */«. ) will be 

solved for directly by taking full advantage of the linearity of 
the solution. The solution starts by observing that the four unknown 
functions (gj^, 8j, fj^, fj) may be written 




- 

1 


0 


0 


0 


3 r 

0 

^ 3 i 

1 

+ fg 

0 

+ f I 

0 

fR 


0 


0 


1 


0 

1 

1 


0 


0 


0 

•m. 


1 


(60) 


Using linearity, the starting values defined by the column matrix 
at X = 0 , 


1 


’ a.. ' 

0 



0 

in+e.gra+'es to the value 

21 

0 


34. 

L J 


Consider first the downstream case. Combining Eqn. (60) with 
the unknown initial conditions described by Eqns. (51), (52) and (55), 
the solution to Eqn. (45) may be written 


1 5 



(\-oC)a^ 

(l - 

(l-ot)a4, 


**“ (l"**®^^*!* + P^I4 
- pa,i + + 

-pa« + 0 + <^)35, + 
-pa42 + (l+ot)a4» + 


= r 


= -a 


= Y 



( 61 ) 


where the quantities a, 3, y and 5 are defined as 
OC = To cos^V, , p = Y = tos^P^j 5 = raSin<f'i (52) 


The solution to the four unknowns a, 3, y» 6 follows inunediately by 
rewriting Eqn. (61) as 


- 


- 



- 

-1 


OC 


a„ ^ 

(a,4 — 3,2.') 

-1 

0 


(311+3,3) 

p 


(aa»-3a.) 

^24 “ ^22) 

0 

1 



Y 


~ ^*1 ) 


-1 

0 


(a^i +333) 

s 


(^4J~ ®4l) 

(344 - ?^42) 

0 

1 


(34, +343) 


This approach is very efficient. It avoids the shooting approach 
used by Davis and Johnson* and the particularly restrictive approach 
used by Eisenberg and Kao\ It is also more accurate that the char- 
acteristic method used by King and Karamcheti®. The work of Eisenberg 
and Kao consists of transforming the coupled acoustic equations 
governing the propagation of sound in a variable area, one dimensional 
duct into two separate equations containing variable coefficients. By 
an appropriate choice of duct shape, the variable coefficients are 
made constant permitting relatively easy solutions of the separate 
equations for particular duct shapes. Physical variables (such as 
pressure and velocity perturbations) are found by tracing back all 
functional transformations. A serious draw back to their approach is 
the severe restriction in duct shape. King and Karamcheti also adopted 
the quasi-one-dimensional assumptions in their solution of the behavior 
of sound in flow. They used the method of characteristics to solve for 
the sound field. Because of the variation of mean flow quantities in 
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the duct, the characteristics are curved. An iterative scheme of 
second order accuracy in mesh spacing is used to account for curva- 
ture when determining the intersection of characteristics of differ- 
ent families from adjacent mesh points. This scheme encounters 
difficulty in maintaining reasonable accuracy. 


To verify that the solutions given by Eqn. (63) are correct, the 
conservation of acoustic energy flux will be independently evaluated 
at X = 0 and 2(KsL). Using the form for the energy flux W given by 
Cantrell and Hart.® 


W(x) = T A(x) 


I II I \f 

(l + M\x))|pV)||u'(x)|cos(0p-0u)+MCx)^-j^^^ 



(64) 


where 0,) and 0.^. represent the phase of the acoustic pressure and 
velocity respectively. 


^ Evaluating Eqn. (64) at x = 0 and x = 2(KgL) and denoting by 
Wi ,Wi-‘ and Wt^ the incident energy flux (at x = 0) the reflected 
energy flux (at x = 0) and the transmitted energy flux (at x = 2(KsL) 
respectively (for the downstream case) , then these quantities may 
be written as 


Wi — 2 . Ao Co I lusl (l + Mo) ^55^ 

Wr “ T Ao Co* 1 |>r I 1 (l - Mo) (66) 

\N^ - ^ AoCoMp(llu^1(l + Mo) (67) 

where = M(x=0) = M(x=2K^L) by symmetry. 

Conservation of acoustic energy requires that 

Wj^ = Wr + (68) 


Upon substitution of the various quantities into Eqn. (68), the 
desired form of the solution is 
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Ref lec.-^eci £ncrau TransmV^eii Enerou 
Coefficvewt (yC^ Coefficie \-t 


where from Ecn. (62), I F,! = oc% p* and 1 Fjl = 

The solution of the upstream case is quite similar to that of 
the downstream case. By combining Eqn. (60) with the boundary condi- 
tions for the upstream case defined by Eqns. (57), (58) and (59), 
the solution is 



Again using the definition of Cantrell and Hart, the conser- 
vation of acoustic energy flux may be written as 



Numerical results for selected duct configurations having throat 
to inlet area ratios (v) equal to 0.305, 0.420, 0.558 and 0.746 and 
for the case of K^L = 1 are summarized in Figures 3-6. The mean flow 
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velocity distribution through the du^L is specified by Eqn. (26). 
Figure 3 suiamarizes the upstream sound propagation solutions and 
Fig. 4, the downstream solutions. The solutions are presented in 
terms of the sound energy reflection coefficient R and the total non- 
dimensionalized acoustic energy flux R + T. In Figs. 3 and 4, energy 
is conserved (i.e., T + R = 1) for throat Mach numbers less than about 
Mth ^ 0.6; this is true for all duct area ratios considered. However, 
for > 0.6, it is clear that acoustic energy flux T + R is not 
conserved, the departure of T + R from unity increasing as the throat 
Mach number approaches unity. 

The departure of T + R from unity for M^-j^ > 0.6 is due to the 
singular behavior of the linearized sound propagation equations. For 
values of near unity, the linearizv ^ sound propagation equations 
are inadequate to account tor the conse ation of acoustic mass and 
momentum. 

Figures 3 and 4 also show that the sound energy reflection 
coefficient (and hence the transmission coefficient) is very sensitive 
to the duct throat to inlet area ratio v. For a fixed value of M^, , 
the increase of R (decrease of T) with the decrease of v is relatea 
to sound reflection or scattering generated by the axial velocity 
gradients of the duct mean flow. This is easily shown by restricting 
the duct mean Mach number at the throat to values less than 0.6 so 
that the effects of compressibility can be ignored. It is straight 
forward to show that the duct mean flow velocity gradient defined 
below as 


V-kW-V,N _ ^ 

L L / 

depends onl' upon the parameter v for a fixed duct geometry and throat 
speed. The large increase in R observed for both the upstream (Fig. 3) 
and downstream (Fig. 4) cases of sound propagation suggests that the 
reduction in sound transmission arises primarily because of scattering 
effects generated by the duct mean flow velocity gradients. This is 
more clearly shown in Fig. 5 which combines the reflection coefficient 
R from Figs. 3 and 4. Figure 5 shows that for values of M^j^ < 0.6, 
there are large and significant differences for the various duct throat- 
to- inlet area ratios shown. The effect, however, on R of upstream or 
downstream sound propagation is small. These results hc'"’^" important 
application to our understanding of the manner in which .j sonic 
inlet acts as an effective barrier to internally generated sound. The 
contribution of the mean flow velocity gradients within the duct may 
contribute more to reduction of sound transmission tiian was previously 
thought . 

Figure 6 is a rather dramatic example of the effect of the duct 
mean flow on the sound field. The distribution along the duct of the 
gradient of tht., imaginary part of the density fj' is plotted for three 
different throat Mach numbers. The singular nature of the sound field 
is evident. 
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4.0 NONLINEAR NUMERICAL EXAMPLE 

The numerical study undertaken in Section 3 showed that the 
solution of the linearized sound propagation equations failed to 
conserve acoustic energy flux (T + R) for throat Mach numbers 
Mth ^ 0.6. The linearized solution is extended below to include 
the previously ignored nonlinear terms. The purpose here is to 
determine numerically, for the test case tieated in Section 3, if 
the accuracy of the solvtion can be extended up to or near = 1 
by including nonlinear terms. 

4 . 1 Derivation of Governing Equations 

The derivation assumes that the mean flow is independent of 
the sound field so that 

= 0 

and 

pViV + ^ =0 

\ dx dx 

With this understanding, the nondimensional sound propagation equa- 
tions are 



(72) 


and 



1 ^ 

ax Q dx 





p 



4 - 


xGi/ S dx 



+ 


u 






( 73 ) 


V. 
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where € represents approximately the ratio of the ampli- 

tude of the incident source sound pressure to the mean flow static 
pressure at x * 0. Equations (72) and (73) have been written with 
the nonlinear terms on the RHS. To solve these equations, we seek 
solutions of the form 



(74) 


and 

u(x . h) = ^ n G )e 

ns-N 


(75) 


The idea behind the assumed form of the solutions is that 
given a sound wave with radian frequency u, nonlinear effects will 
generated higher and lower harmonics of frequencies ranging from 
-No), - (N-l)a:, . . . -2u», -u, o, o>, 2u,... (N-l)u, Nu) for fixed N. Since 
we are interested only in determining the importance of retaining 
the nonlinear terms, we select N = 2. Thus we will have to solve 
five simultaneous differential eqns. with variable coefficients. 


To demonstrate the details of the proposed approach, we will 
consider the first terms on the LHS and RHS of Eqn. (72). 




Substituting Eqns. (74) and (75) yields 

2 2 2 . 


'(•^0 ^ n. 


in't 


n*-2 

Multiplying both sides by e 
2 -T 


(j*K)t 


K»-2 j=-2 


-imt 






2 2 


and integrating over a period T yields 

T 

(76) 


K »-2 




r 

e 


All terms on the LHS vanish except the tern n - m and similarily all 
combinations of terms on the RHS vanish ::xcep those summing to m 
(i.e., 3 + k « m) . Thus Eqn. (76) may je written 
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To compare directly with the linear case, we select m = 1 so that the 
contribution of the nonlinear term to the fundamental solution is 



Thus to correct the linear solution for nonlinear effects requires 

the solution to the functions f and for all n up to +2. 

n °n — 

The functions fo and go require special comment. They correspond 
to time independent or steady state solutions (the so-called acousti- 
cal steaming solutions). We believe that these solutions require 
at least two -dimensionality to exist because the streamlines are 
closed (i.e., there is no mass addition). Thus they vanish in the 
present one-dimensional application. 

Incorporating the nonlinear terms, the governing differential 
equation becomes 




— - 

— - 


-a 

IVl A A T|r,n 

3',n 

_ $ Cn-0 

M. 


6 

1 


^r.n 


L 


l-M* 

^ ^ ^r,n 

Kr. 



■fi.n 



_ 


_ 



_ * 


(79) 


where the matrix [MlI is the same as the linear matrix defined by 
Eqn. (45), 6(n-l) is the kroncker delta, and the quantities S and T 
are nonlinear source terms defined as 







( 80 ) 






Z E 

iw K 




^ m 



± f 

■*■ R dx VR/ 3'**'“ 


( 81 ) 


Subscripts m and k denote respectively the m-th and k-th mode of the 
sound wave. The subscript n ranges from -2 to +2. It is important 
to realize that the solution to Eqn. (79) requires that the functions 
fn and gn be known for all n (in the present application n = *2). 

4 . 2 Boundary Conditions 

The boundary conditions specified in Section 3.2 for the linear- 
ized sound propagation equations can be extended, in a straight-forward 
manner, to apply to the nonlinear equation described by Eqn. (79). 

The boundary conditions are different for the cases of downstream and 
upstream sound propagation. 


4.2.1 Downstream Propagation 


At x=0, both incident and reflected waves exist. The incident 
wave is assumed known, of frequency corresponding to n = 1, and 
generated far upstream of the constricted. Reflected waves having 
(complex) frequencies corresponding to n = -2, -1, 2 are generated by 
the constriction. Also, transmitted waves are generated at x = 2KgL 
corresponding to values of n = -2, -1, 1, 2. The corresponding func- 
tions fn(0) and gn(0) are 


and 



S(n-0 


-H 


O.rt 




Incident Wave 


Reflected Wave 


(82) 



V 


S(n-l) ^ 


V-. 


fo.n 



(83) 


Incident Wave 


Reflected Wave 


where 6(n-l) is used to specify the non-dimensional source sound in- 
troduced far upstream of the duct constricted. The corresponding 
functions f„ and ^„(2KsL)are 

= 3"(2Ks 1-) = e''*’*’’' (84) 


23 



where 




Z n KsL 
l + M„ 


C85) 


The equation describing the conservation of energy flux for 
the downstream propagation case is, 

H 


I 

n=-N 


. V 3 3_ 



ir 1 1 

\ H-M. / ' ' 

« •‘•a.n • 

> V , / V 

y ■/ 


= 1 


( 86 ) 


Reflected 
4.2.2 Upstream Propagation 


Transmitted 


Analogous to the downstream sound propagation, the upstream 
boundary conditions are 

2i 

= S(n-l)e + 

V / 






(87) 


Incident Wave 


Reflected Wave 


and 




^ / N A* 

~d(«-l)e l-Mo 






2.n 


( 88 ) 


Incident Wave 


Reflected Wave 


where 




2n L 

l-M. 


(89) 


The equation describing the conservation of energy flux, for the 
upstream sound propagation case, is 
N 


r ■ ir.,:* 

Z ' =-y ' ' ' 




(90) 


4.3 


n=-N .Reflected 

Integration Scheme 


Transmitted 


The value of N is restricted to N = 2 . Consider first the 
downstream case. There are a total of 16 unknowns 
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corresponding to values of N of (-2, -1, 1, 2). From Eqn. (79), 
there are 16 equations to be solved. To solve these equations, the 
nonlinear source terms are assumed known. Thus the governing 
equations are linear and the general solution is a linear combination 
of the homogeneous and particular solutions for each value of n. The 
solution proceeds much like that described in Section 3 except that the 
contribution to the particular solution must be included. For the 
downstream case, the n^h solution may be written 



(91) 


where 








“ *^),n ~ (92) 


and 

. 3r% • f M represent the particular solutions of g and f. 
The solution to Eqn. (91) for values of n = -2, -1, 1, 2 yields the 
values of the 16 unknowns. The solution may be written as 


f- 


■ 



— 

-1 

p 




W-as) 

-1 

o 
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The solution to the upstream sound propagation case is similar 
to that derived above and is written below as 


25 



r *1 




m 

-1 

fa 



K-»") 
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0 1 





(a; -a.") 

(a»-a,I) 

-1 0 





-a*) 

K-a«) 

0 I 




(94) 


In order for this scheme to be successful, the various source 
terms for Sn and Tn must be known for each value of n. The following 
approach is used to solve for Sn and Tn: 

(1) Initially, the homogeneous solution n = 1 is solved. From Eqns. 
(80) and (81) , the source terms for n = 2 is seen to depend only upon 
the n = 1 solutions. This follows because fo and are identically 
zero leaving only the terms m = k = 1 summing to m = k = 2) . Thus 
the n = 1 homogeneous solution is used to estimate the f 2 and g 2 
solutions . 

(2) With the homogeneous n = 1 and the nonhomogeneous solutions known, 
the n = -2 and n = -1 solutions remain to be estimated. The nonlinear 

n = -2 solution is estimated directly from the nonlinear n = 2 solution 
by assuming that Re(fjj) = Re(f.jj) and Im(fn) = Im(f-n) with similar 
expressions for gn . This assumption follows from the physical con- 
straint that the solutions for f^ and gn must be real. Thus knowing 
the solution for n = +2, the solution for n = -2 is also known. 

(3) The solution for n = -1 can be estimated from the n = -2 and 

n = 1 solutions. This follows from Eqn. (79) where now m + k = -1; 
the only combination permitted is m = -2 and j * 1. 

(4) With the solutions for the cases n = -1, 2 and -2 known, the 
correction to the linear (i.e., homogeneous) n = 1 solution is obtained 
from Eqn. (79) where here the various combinations summing to m + j = 1 
are m = 2, j = -1, and m = -1, j = 2. 

(5) Steps (1) through (4) above should be repeated until satisfactory 
convergence is obtained. 

4.4 Results 


The results of this approach are summarized in Figures 7 and 8 
for the cases k^L = 1 and v = 0.746. Figure 7 represents the down- 
stream sound propagation solution and Figure 8, the upstream solution. 
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The results are discouraging in the sense that none of the solutions 
adequately conserved acoustic energy flux. 

Figures 7 and 8 represents the results of only a single integra- 
tion pass. No iterations were attempted because of contract cost 
limitations. Thus convergence was not explored, the accuracy of the 
proposed solution has not been adequately pursued. For this reason, 
there will be no further interpretation of the nonlinear solution. 

This does not mean that the approach is invalid -- it only means that 
more than one iteration is required before satisfactory convergence 
is achieved. 

5.0 CONCLUDING REMARKS 

The linearized equations describing the propagation of sound 
in variable area ducts have been shown to be singular when the duct 
mean axial flow component is sonic. The singularity is removed when 
previously ignored nonlinear terms are retained. This conclusion is 
quite general because the flow field within the duct is fully three- 
dimensional (both mean and fluctuating parts) and unsteady. The only 
assumptions made are that the fluid is inviscid and non-heat conduct- 
ing. 


A numerical study was conducted to map out, for the case of a 
plane wave propagating in a one-dimensional converging-diverging duct, 
the axial mean flow Mach number range at the throat M-^}^ for which the 
linearized sound propagation equations are valid. The results of the 
study showed, in terms of the conservation of acoustic energy flux, 
that the linearized sound propagation equations are valid only when 
the mean flow Mach number M^jj < 0.6 at the duct throat. For values of 
Mth > 0.6, the acoustic energy flux is not conserved. Since consider- 
able care has been taken to minimize numerical errors by varying in- 
tegration step size, integration algorithm, etc., it is concluded 
(at least for this example) that for M^)^ > 0.6, the linearized sound 
propagation equations do not adequately model the conservation of 
acoustic mass and momentum within the duct. 

The numerical study also suggested for Mth 0.6, that sound 
transmission within a converging-diverging duct is reduced primarily 
by scattering due to the duct mean velocity gradients. This does not 
preclude the possibility that other velocity gradients, say in the 
radial direction, may also significantly reduce sound transmission. 
Recall that this study considers only plane waves propagating in a 
one-dimensional mean flow. 

These conclusions have an important bearing on our present under- 
standing of how the sonic inlet acts as an effective noise barrier to 
internally generated sound. The large observed decreases of the sound 
transmission in a sonic inlet may arise in large mec.-ure because of 
scattering due to mean flow gradients within the inlet. Previously, 
it was believed that most of the sound transmission reduction was due 
to the mean sonic flow at the throat opposing, thereby preventing, 
the propagation of sound through the throat. 
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An attempt was made to study the nonlinear behavior of sound 
waves propagating in variable area ducts containing sonic or near 
sonic flow. Meaningful results were not obtained. This does not 
necessarily mean that our approach is invalid, but rather that more 
effort is required to assess its validity. It would be of immense 
value to continue the nonlinear solution described herein. This 
would provide a preliminary understanding of the acoustic behavior 
of the sonic inlet. 
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FIGURE 2(a) DOWNSTREAM SOUND PROPAGATION 



FIGURE 2(b) UPSTREAM SOUND PROPAGATION 
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